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1.  Introduction 


The  purpose  of  this  report  is  to  set  forth  a  new  set  of  cirrus  cloud  models  for  use  with  the  MOD- 
TRAN4  atmospheric  transmission  model  and  to  provide  a  wider  range  of  cirrus  models  with  which  to 
investigate  future  remote  sensing  techniques  that  may  have  to  observe  through  cirrus.  The  cirrus 
cloud  models  were  derived  by  selecting  mode  particle  sizes  between  and  beyond  the  two  existing  cir¬ 
rus  cloud  models  in  MODTRAN4:  Standard  (mode  particle  radius  64  pm)  and  subvisual  (mode  4 
pm.)  The  new  set  of  models  consists  of  70  models  whose  mode  particle  radius  range  from  between 
0.1  and  64  pm.  The  format  of  the  models  meets  the  input  requirements  for  MODTRAN4  and  consist 
of  files  containing  the  input  information  on  card  2E2.  All  the  model  calculations  are  based  on  spheri¬ 
cal  particles.  A  plan  is  included  for  extending  the  models  to  more  realistic  ones  that  use  nonspherical 
particles. 

The  Moderate  Resolution  Transmittance  (MODTRAN)  code  developed  by  the  Air  Force  Research  ^ 
Lab  (AFRL)  calculates  atmospheric  transmittance  and  radiance  for  frequencies  from  0  to  50,000  cm 
at  moderate  spectral  resolution,  primarily  2  cm  (20  cm  in  the  UV)  (Berk  et  al.  1989).  The  devel¬ 
opment  of  the  MODTRAN  model  was  motivated  by  the  need  for  higher  spectral  resolution  than  was 
available  in  the  Low  Resolution  Transmittance  (LOWTRAN7)  (Kneizys  et  al.  1986).  MODTRAN’s 
capabilities  include  spherical  refractive  geometry,  solar  and  lunar  source  functions,  scattering 
(Rayleigh,  Mie,  single  and  multiple),  and  default  atmosphere  profiles  (gases,  aerosols,  clouds,  fogs, 
and  rain).  MODTRAN,  version  4,  release  1  was  used  for  these  calculations,  and  is  the  most  current 
release. 

The  Cirrus  cloud  models  in  MODTRAN,  and  its  predecessor  LOWTRAN,  were  calculated  by  Shettle 
et  al.  (1988)  using  the  optical  constants  for  water  ice  measured  by  Warren  (1984).  These  are  known 
as  the  “Standard”  and  “Subvisual”  Cirrus  cloud  models.  The  shape  of  the  ice  crystals  was  approxi¬ 
mated  using  spherical  particles,  allowing  the  use  of  Mie  theory  to  calculate  the  absorption  and  extinc¬ 
tion  (Bohren  1983,  van  de  Hulst  1957).  Shettle  et  al.  assumed  a  particle  size  distribution  given  by  the 
log-normal  distribution  function: 


dN/(Lr  =  n(4)  =  aaexp(-ba),  ( 1 ) 

where  ais  the  particle  radius,  and  a  and  b  are  constants.  Figure  1  shows  the  particle  size  distributions 
used  to  compute  Kg  for  the  standard  and  subvisual  models.  The  mean  (mode)  particle  radius  for  the 
standard  and  subvisual  models  is  64  |xm  and  4  pm,  respectively,  and  the  calculations  were  done  under 
the  assumption  that  all  the  particles  are  spherical.  The  particles  were  treated  individually  between 
0.05  pm  and  1000  pm  with  799  interspaced  sizes. 

The  term  “cirrus”  is  used  to  describe  a  number  of  clouds  with  similar  characteristics  (Lynch  and  Sas- 
sen  2001),  the  main  difference  within  the  cirrus  category  being  particle  size.  In  a  previous  study. 
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Particle  Size  Distributions  for  Cirrus 


Radius  a  Oim) 

Figure  1 .  The  particle  size  distribution  for  the  subvisual  model  contains  more  particles 
cm  ,  and  they  are  on  average  much  smaller  (^mode  =  4  pm)  than  for  the  stan¬ 
dard  model  ((^rnode  =  64  pm). 

Lynch  and  Mazuk  (2000)  discussed  how  the  particle  size  distribution  has  a  strong  influence  on  the 
absorption  and  emission  spectra  for  a  particle  cloud.  The  paper  by  Fu  (1996)  summarizes  25  distinct 
particle  size  distributions  that  have  been  experimentally  measured  in  cirrus  clouds.  Fu  also  notes  that 
the  particle  size  distributions  are  difficult  to  measure  experimentally.  Small  ice  particles  are  particu¬ 
larly  difficult  to  measure,  and  in  some  experiments,  the  distribution  functions  were  extrapolated  to 
small  sizes  from  measurements  of  larger  particles. 
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2.  Derivation  of  the  New  Models 


The  new  cirrus  models  fill  in  the  mode  particle  sizes  between  the  two  existing  cirrus  models  and 
extend  the  subvisual  cirrus  model  to  smaller  particle  sizes  than  have  previously  been  available.  The 
same  modeling  methods  used  by  Shettle  et  al.  (1988)  were  used  here.  The  Mie  scattering  code  of 
Shettle  et  al.  was  used  along  with  the  optical  constants  of  Warren  (1984)  to  generate  a  set  of  extinc¬ 
tion  coefficient  (Ke,  km"\  absorption  coefficient  (Ka,  km"  ),  and  asymmetry  parameter  (g),  all  as  a 
function  of  wavelength  between  0.4  and  14.7  pm. 
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3.  The  New  Models 


The  new  cirras  cloud  models  consist  of  a  large  amount  of  numerical  tabulations  of  extinction, 
absorption,  and  asymmetry  coefficients.  Therefore,  we  show  them  here  in  graphical  form  in  Figures 
2-15.  The  actual  numerical  models  will  be  in  the  form  of  MODTRAN  card  2E2  (See  below)  and  will 
be  transferred  electronically.  In  each  of  Figures  2-9,  the  standard  cirrus  model  corresponds  to  the  far 
edge  of  the  plot  (%ioDE  =  64  pm),  and  the  subvisual  cirrus  model  corresponds  to  the  near  edge  {a^  = 
4  pm).  In  each  of  Figures  10-17,  the  subvisual  cirrus  model  corresponds  to  the  far  edge  of  the  plot 
(4m  =  4  pm). 

Figure  2  shows  a  surface  plot  of  the  extinction  coefficient  Kg  over  the  wavelength  range  of  0.4  to 
14.7  pm.  This  is  normalized  to  unity  at  0.55  pm  for  all  mode  particle  sizes,  but  can  be  rescaled 
according  to  the  user’s  selection  of  optical  depth.  Note  that  Kg  is  relatively  flat  and  unstructured 
except  in  the  foreground  where  the  small  particle  effects  become  evident.  The  two  obvious  valleys 
correspond  to  the  vibrational  structure  of  ice  at  3.16  pm  and  the  broader  rotational  structure  around 
10.5  pm. 

Figure  3  shows  a  surface  plot  of  the  absorption  coefficient  Ka  over  the  wavelength  range  of  0.5  to  14 
pm.  Here,  there  is  considerably  more  stmcture  in  both  wavelength  and  mode  particle  size  space,  and 
as  in  Figure  2,  the  structure  increases  with  decreasing  particle  size. 

Figure  2  is  relatively  flat  because  the  extinction  is  dominated  by  scattering,  and  scattering  efficiency 
is  not  strongly  dependent  on  particle  size  or  wavelength.  In  fact,  it  depends  on  the  scattering  size 
parameter  X  =  litalX,  which  is  generally  greater  than  unity  for  most  of  the  situations  shown  in  the 
figure.  The  ample  structure  in  Figure  3,  on  the  other  hand,  is  the  result  of  numerous  variations  in  the 
optical  constants  of  ice.  These  are  embodied  in  the  absorption  size  parameter  Q  =  47ck4' /X  (Lynch 
and  Mazuk  1998),  where  k  is  the  imaginary  part  of  the  index  of  refraction.  In  this  situation,  Q,  «  1, 
and,  therefore,  the  variations  in  k  become  dominant. 

The  asymmetry  parameter  g  results  are  shown  in  Figure  4,  showing  similar  structure.  There  is  con¬ 
siderable  structure,  but  unlike  the  values  in  Figure  3,  g’s  dynamic  range  is  fairly  small,  generally  fal¬ 
ling  in  the  range  of  0.8  to  1 .0.  As  4m  decreases,  the  particles  become  small  compared  to  the  wave¬ 
length,  and  the  X  drops  below  unity.  When  this  happens,  the  particles  become  more  like  “Rayleigh” 
particles,  and  g  approaches  zero. 

Figures  5,  6,  and  7  are  expanded  version  of  Figure  2,  3,  and  4.  These  figures  show  the  extinction, 
absorption,  and  asymmetry  parameters  between  1  and  1 .5  pm.  In  Figure  5,  the  extinction.  Kg,  is 
relatively  flat,  although  some  size-  and  wavelength-dependent  “warpage”  is  present.  Figure  6  shows 
virtually  no  variation  in  the  absorption  Ka,  but  Figure  7  shows  significant  variation  in  g  as  a  function 
of  mode  radius  4m- 
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Figures  8  and  9  show  the  scattering  coefficient  Ks,  which  is  defined  as  Kg  -  Ka-  They  are  included 
for  completeness. 

It  is  clear  from  all  of  the  above  results  that  changes  in  the  mean  particle  size  (i.e.,  changes  in  ) 
have  the  greatest  influence  for  small  particles.  This  sensitivity  is  well  known  and  has  been  discussed 
in  the  context  of  cirrus  clouds  by  Takano  et  al.  (1992).  Therefore,  we  extended  the  range  of  values 
for  downward  to  between  0.1  and  4.0  pm  and  calculated  Kg,  Ka  and  g.  The  results  are  shown  in 
Figures  10-17,  and  correspond  to  Figures  2-9,  except  they  are  for  smaller  particles.  In  each  of  Fig¬ 
ures  10-17,  the  subvisual  cirrus  model  corresponds  to  the  far  edge  of  the  plot  (a^  =  4  pm). 


Figure  8.  Scattering  coefficient  (Ke  km  '):  X  -  0.5 
to  14  pm,  and  Om  =  4  to  64  pm. 
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Figure  9.  Scattering  coefficient  (Ke  km  ^):  X  =  1.0  to 
1 .5  [im,  and  ay[  =  4  to  64  ^im. 


Figure  10.  Extinction  coefficient  (Ke  km  ^):  X  =  0.5 
to  14  [im,  and  a'M  =  0A  to  4  [xm. 
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4.  Using  the  New  Models  in  MODTRAN4 


The  new  cirrus  cloud  models  were  developed  for  use  with  the  MODTRAN4  atmospheric  transmission 
model.  MODTRAN4  allows  for  the  insertion  of  a  cloud  model,  which  is  then  used  in  conjunction 
with  model  (or  measured)  atmospheric  data  in  the  calculation  of  transmission  and  radiance.  Inputs  to 
MODTRAN4  consist  of  lines  of  formatted  ASCII  text,  and  are  referred  to  as  “cards”  in  reference  to 
IBM  punch  cards  formerly  used  as  computer  inputs. 

The  cloud  and  aerosol  specification  begins  on  Card  2.  Here,  we  commonly  set  the  parameter  IHAZE 
=  1  for  a  rural  aerosol  model  with  23  km  visibility.  We  also  set  ICLD  =  1  to  select  a  cumulus  cloud. 
Remaining  values  are  set  to  zero,  but  the  user  can  make  further  selections  depending  on  the  charac¬ 
teristics  of  a  particular  analysis. 

Card  2A  (alternate  form  used  if  ICLD  =  I  through  10)  defines  the  parameters  for  the  cloud  model 
chosen  on  Card  2,  and  is  used  for  clouds  other  than  cirrus.  Entries  on  this  card  trigger  the  reading  of 
cards  2E1  and  2E2.  The  simplest  way  to  define  a  cirrus  cloud  is  to  specify  its  base  height  and  thick¬ 
ness  (both  in  km)  and  its  extinction  coefficient.  Kg,  from  the  calculated  models.  In  this  mode,  how¬ 
ever,  the  thickness  must  be  chosen  so  as  to  give  the  proper  optical  depth  at  0.55  [xm.  The  optical 
depth,  T,  and  Kg  are  simply  related  through  the  cloud  thickness. 


x  =  KeL, 

where  L  is  the  cloud  thickness  in  km.  The  user  also  sets  the  number  of  layer  boundaries  used  to 
define  the  cloud  structure  on  card  2E1,  and  the  number  of  wavelengths  to  be  input  for  cloud  extinc¬ 
tion  and  absorption  values  on  card  2E2.  The  cloud  column  water  density  is  set  to  zero.  Finally,  the 
user  specifies  the  model  asymmetry  values  rather  than  defaulting  to  Henyey-Greenstein  values. 

3 

Card  2E1  defines  the  cloud  layer’s  ice  and  liquid  water  densities  in  g/m  .  As  it  happens,  the  value  of 
densities  are  not  used  directly;  so  as  long  as  they  are  greater  than  zero,  the  code  will  work.  For  cirrus 
clouds  only,  the  user  must  zero  out  the  water  density  and  the  rain  rate  as  well. 

Card  2E2  defines  the  extinction,  absorption,  and  asymmetry  of  water  and  ice  for  each  wavelength 
(Table  1).  Here,  the  user  must  input  negative  values  for  the  water,  which  causes  the  cloud  model  to 
use  the  default  extinction,  absorption,  and  asymmetry  values  for  the  cloud  we  specified  in  Card  2 
(here  a  cumulus  cloud).  But  recall  that  the  user  previously  zeroed  out  the  water  content  of  the  cloud, 
so  there  is  no  contribution  here.  So  the  only  values  that  should  be  active  are  the  ones  input  for  the  ice 
based  on  the  models  presented  here. 
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Table  1 .  Sample  input  Card  2E2 


Units 

1.38 

Wavelength 

(pm) 

-9 

Extinction  coefficient 

(km  ^ 

for  water  (negative  suppresses  water  cloud) 

-9 

Absorption  coefficient 

(km  ) 

for  water  (negative  suppresses  water  cloud) 

-9 

Asymmetry  parameter 

(unitless) 

(negative  suppresses  water  cloud) 

0.1 

Extinction  coefficient 

(km  ^ 

for  ice  (typical) 

0.01 

Absorption  coefficient 

(km  ^ 

for  ice  (typical) 

0.86 

Asymmetry  parameter 

(unitless) 

for  Ice  (typical) 
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5.  Limitations  of  the  Models 


The  models  presented  here  are  based  on  Shettle  et  al.’s  cirrus  cloud  models.  They  give  reasonable 
results  and  reproduce  the  basic  properties  of  cirrus  clouds,  i.e.,  absorption  peaks  at  the  appropriate 
places,  representative  transmission,  scattering  cross  sections  that  behave  as  expected,  etc.  To  the 
extent  to  which  the  two  models  in  MODTRAN4  are  useful,  the  models  presented  here  extend  this 
utility  to  intermediate  size  regimes  and,  perhaps  more  importantly,  to  smaller  sizes  that  are  associated 
with  subvisual  cirrus.  For  most  applications,  they  are  about  as  good  as  can  be  achieved  within  the 
computational  limitations  of  MODTRAN. 

All  the  new  cirrus  cloud  models,  as  well  as  the  built-in  cirrus  cloud  models,  suffer  from  one  major 
limitation:  they  are  calculated  assuming  that  the  clouds  are  composed  of  spherical  particles.  The  ice 
particles  comprising  cirrus  clouds  are  unquestionably  non-spherical,  often  highly  so.  The  impact  of 
the  spherical  particle  assumption  is  that  the  scattering  of  light  is  not  accurately  calculated  over  all 
angles;  i.e.,  the  bi-directional  reflectance  distribution  function  (BRDPO  is  not  well  modeled.  How¬ 
ever,  the  current  release  of  the  MODTRAN4  model  does  not  allow  insertion  of  a  BRDF  for  any  cloud 
model;  only  ground  surface  BRDFs  can  be  inserted,  so  this  solution  is  not  available.  The  alternative 
is  to  parameterize  the  effects  of  non-spherical  particles  within  the  limitations  of  the  MODTRAN 
inputs.  This  approach  will  be  explored  in  the  next  phase  of  our  modeling  efforts. 
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6.  Conclusions 


We  have  produced  a  new  set  of  cirrus  cloud  models  for  use  with  the  MODTRAN  atmospheric  trans¬ 
mission  code.  These  models  extend  the  range  of  mean  particle  sizes  beyond  the  limits  of  the  existing 
models  included  in  MODTRAN,  covering  the  broader  range  of  cloud  composition  seen  in  nature. 
These  models  can  be  used  in  conjunction  with  MODTRAN  to  investigate  the  effects  of  cirrus  clouds 
on  future  remote  sensing  techniques  that  may  have  to  observe  through  cirrus. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  spe¬ 
cializing  in  advanced  military  space  systems.  The  Corporation’s  Laboratory  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staff’s  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analy¬ 
sis,  solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and 
CCD  detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid 
state  laser  design,  micro-optics,  optical  communications,  and  fiber  optic  sensors;  atomic 
frequency  standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation 
and  beam  control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evalua¬ 
tion,  battery  electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  fail¬ 
ure  analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aero  thermodynamics;  chemical  and 
electric  propulsion;  environmental  chemistry;  combustion  processes;  space  environment 
effects  on  materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and 
structural  control;  lubrication  and  surface  phenomena. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmos¬ 
pheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  sur¬ 
veillance,  imaging,  remote  sensing,  and  hyperspectral  imaging;  effects  of  solar  activity, 
magnetic  storms  and  nuclear  explosions  on  the  Earth's  atmosphere,  ionosphere  and  magne¬ 
tosphere;  effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation,  design  fabrication  and  test;  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific  chemical 
reactions  and  radiative  signatures  of  missile  plumes. 

Center  for  Microtechnology:  Microelectromechanical  systems  (MEMS)  for  space 
applications;  assessment  of  microtechnology  space  applications;  laser  micromachining; 
laser-surface  physical  and  chemical  interactions;  micropropulsion;  micro-  and  nanosatel¬ 
lite  mission  analysis;  intelligent  microinstruments  for  monitoring  space  and  launch  sys¬ 
tem  environments. 

Office  of  Spectral  Applications:  Multispectral  and  hyperspectral  sensor  development; 
data  analysis  and  algorithm  development;  applications  of  multispectral  and  hyperspectral 
imagery  to  defense,  civil  space,  commercial,  and  environmental  missions. 


